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Qh. Abstract 

Nonequilibrium dynamics in quantum field theory has been studied exten- 
sively using truncations of the 2PI effective action. Both 1/N and loop ex- 
pansions beyond leading order show remarkable improvement when compared 
5_i ■ to mean-field approximations. However, in truncations used so far, only the 



leading-order parts of the self energy responsible for memory loss, damping and 
equilibration are included, which makes it difficult to discuss convergence sys- 
tematically. For that reason we derive the real and causal evolution equations 
for an O(N) model to next-to-next-to-leading order in the 2PI-l/iV expansion. 
Due to the appearance of internal vertices the resulting equations appear in- 
tractable for a full-fledged 3 + 1 dimensional field theory. Instead, we solve 
the closely related three-loop approximation in the auxiliary-field formalism 
numerically in + 1 dimensions (quantum mechanics) and compare to previous 
approximations and the exact numerical solution of the Schrodinger equation. 
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1 Introduction 

The two-particle irreducible (2PI) effective action formalism has proven very powerful 
for out-of-equilibrium quantum field theory over a wide range of applications pQ. Since 
it necessarily employs an expansion and truncation of the effective action, one should 
be concerned with how well these expansions converge. 1 

An extensive study has been made of effective mean-field approximations (Gaus- 
sian, leading order (LO) 1/N or Hartree approximations) which amounts to including 
a self-consistent, time-dependent effective mass in the dynamics of the one- and two- 
point functions [3]. However, since such an approximation does not include scatter- 
ing, phenomena such as damping, memory loss and equilibration cannot be described 
properly This can be traced back to the existence of (infinitely many) conserved 
charges in the mean-field dynamics and the presence of a nonthermal fixed point 
01 Ej- 2 Going beyond the mean-field approximation, by either employing the weak- 
coupling or the 1/N expansion to next-to-leading order (NLO), effective memory 
loss, damping and equilibration are present. 3 In particular it has been found that 
at late times the system evolves towards a quantum equilibrium state, characterized 
by suitably defined field occupation numbers approaching the familiar Bose-Einstein 
or Fermi-Dirac distribution functions. The self-consistently determined propagators 
have a time-dependent mass and width. It is found therefore that NLO approxima- 
tions improve dramatically upon mean-field approximations and qualitatively repro- 
duce the dynamics expected on physical grounds from the full, untruncated system. 

The natural question to ask is whether a truncation at NLO also gives quanti- 
tatively correct results. This issue has been investigated in several ways, usually by 
comparing the NLO results to other available approximations. For example, when the 
coupling is small, normal perturbation theory should be applicable for e.g. estimates 
of damping rates. Indeed, the results in ref. [12] suggest that the perturbative result 
is reproduced within a factor of two for small coupling, where the difference is due 
to the effect of including self-consistent infinite resummations in the 2PI approach. 
However, this amounts to a test of perturbation theory rather than a verification 
of the 2PI formalism. Similarly, it is well-known [22J how to derive on-shell kinetic 
(Boltzmann) equations from truncations of the 2PI effective action and one can com- 
pare the self-consistent 2PI dynamics with dynamics from kinetic theory [XT ] I23 | I24| . 

1 Of course, perturbative or 1/N expansions in quantum mechanics and field theory are usually 
not convergent but instead asymptotic, see e.g. ref. [2j for discussions concerning the O(N) model. 

2 Some scattering is present if one allows for inhomogeneous mean fields [S1[7||H]. 

3 See e.g. refs. E3 EH U2 for the loop expansion, refs. Ql El HU EU 113 Hi E3 HU for the 
2PI-1/JV expansion to NLO, and ref. J2J for 2PI dynamics with fermions. 



Figure 1: Weak-coupling expansion of the 2PI effective action in the broken phase 
of a scalar \<f) theory ^2] • 

However, this again serves more as a test of kinetic theory than of the 2PI truncation. 
Another possibility, relevant for the dynamics at very late times, is to study trans- 
port coefficients from the 2PI effective action, which gives insight into what scattering 
processes are included |2o"ll2l)] . 

When the exact, untruncated evolution is accessible, one may carry out a direct 
comparison. This option is available in quantum mechanics, where the dynamics 
from the 2PI effective action can be benchmarked against the numerical solution of 
the Schrodinger equation j^Zj. Finally, perhaps the most detailed comparison has 
been made within classical statistical field theory, where direct numerical simulations 
are straightforward jH] . Moreover, the 2PI formalism can be easily applied to classical 
dynamics of an initial nonequilibrium ensemble JTSJE3 120] - 4 I n ref- ^H] it was found 
in the context of the 1 + 1 dimensional O(N) model that the nonperturbative classical 
evolution, obtained numerically, is well reproduced by the 2PI-1/N expansion at NLO 
for iV larger than about 10, providing direct support for the use of the 2PI-1/N 
approximation at NLO. 

Ideally, within a given approximation scheme, the validity of a truncation is es- 
timated by extending the method to the next order in the expansion and test for 
effective convergence. Due to the absence of scattering in mean-field approximations, 
it is essential to go to NLO and compare NLO with approximations that go further. 
In the coupling expansion, a first attempt to do this was made in ref. |T2j. Here 
it was noted that in the broken phase of a A0 4 theory two diagrams contribute at 
O(X): the Hartree and the background field dependent diagram, shown in fig. [T] on 
the left. Perturbatively neither of these diagrams leads to (on-shell) damping, but 
after the self-consistent 2PI resummation the second one does. This approximation 
can then be extended by inclusion of the basketball diagram, shown in fig. Q on the 
right. From a comparison between the dynamical evolution in the two truncations, it 
was found that equilibration times differ, but only by a factor less than two, see ref. 
for more details. Both truncations are at the level of complexity of the 2PI-l/iV 



4 In quantum field theory, stochastic quantization techniques have recently yielded the first real- 
time nonequilibrium lattice simulation results [28]. When developed further, this would offer the 
possibility for direct tests as well. 



expansion to NLO and are therefore numerically tractable. Technically, the common 
feature of these truncations, and in fact of all truncations treated so far, is the absence 
of internal vertices in self energies. This is important from a numerical point of view, 
which is complicated due to the presence of memory integrals in the evolution equa- 
tions. Additional vertices require extra memory integrals to be carried out, which is 
numerically expensive. 

In this paper we extend the analysis and include for the first time self energies with 
internal vertices, using the framework of the 2PI-l/iV expansion to next-to-next-to- 
leading order (N 2 LO) in the symmetric phase of the O(N) model. This opens up the 
possibility to make a quantitative comparison between evolution at NLO and N 2 LO 
in the 1/iV expansion. The N 2 LO effective equations of motion contain additional 
(nested) space-time integrals for each time step, when compared to NLO. It is beyond 
the scope of this paper to solve those equations in full field theory. For the sake of 
illustration we instead specialise to quantum mechanics which allows us to test the 
consistency of the equations, the conservation of energy and effective convergence of 
the expansion. We stress that the dynamics of quantum mechanics is of course very 
different from field theory. Given sufficient computer power the evolution equations 
can readily be implemented in field theory. 

The paper is organized as follows. In the next section we briefly review the 2PI- 
1/N expansion in the O(N) model, following closely the discussion and notation of 
ref. J7|. In section |3] we present the dynamical equations to N 2 LO and give explicitly 
(part of) the statistical and spectral self energies, which are much more involved 
than at NLO due to the presence of internal vertices. Results from the numerical 
implementation for the + 1 dimensional case are shown in section 01 while the 
outlook is given in section In three appendices we collect technicalities related to 
the standard loop expansion, multi-loop contour integrals and the numerical solution 
of the Schrodinger equation. 

2 2PI-1/7V expansion 

Throughout the paper we consider an N component scalar field in the O(N) sym- 
metric phase ((4> a ) = 0). The action is 



S = — J d A x 



-d^a&fa + -m 2 <P a <P a + —{<p a <p a ) 2 



(2.1) 



with a = 1, . . . , N. Doubled indices are summed over. As appropriate to an out- 
of-equilibrium treatment, the fields are defined along the Keldysh contour C in the 



complex-time plane, see appendix [BJ The 2PI effective action depends on the full 
two-point function G of the theory and can be parametrised as J2H] 



1 - ■ ■ -Y-l , l rr„ n-\ 



T[G] = -Tr lnG" 1 + -Tr G, 1 (G - G ) + T 2 [G], (2.2) 

where Gq 1 denotes the free inverse propagator. Variation of the effective action with 
respect to G results in the Schwinger-Dyson equation for the two-point function, 
G^ 1 = Gq 1 — E, which after multiplication with G reads 

- [D x + m 2 ] G ab {x, y) =i S ac (x, z)G cb {z, y) + i5 ab 5 c {x - y). (2.3) 

J z 

Here we used the short-hand notation 

r d 4 z. (2.4) 



'C 

The self energy is given by 

S„, = 2i^M. (2.5) 

oG ab 

The full four-point function is represented by the nonlocal term on the RHS of eq. 
(J2.3J) . For future reference we note here that 

TTjy {4> 2 a(x)(f>l{x)) =i d 4 z S ac (x, z)G ca (z, x), (2.6) 

which can be verified using e.g. the Heisenberg equations of motion. 

To implement the 1/N expansion efficiently, it is convenient to use the auxiliary- 
field formalism J2H1 EH EH EZJ • The action then reads 



S[4>,x} = - 



1 1 3iV 1 

-d^ a d^(l) a + -m 2 <i) a <i) a - —x 2 + - 



(2.7) 



Integrating out \ yields the original action (|2.1|) . The 2PI effective action is now 
written in terms of the one-point function \ = (%(a?)) and the two-point functions 

G ab (x,y) = (TcMx)My)), D(x,y) = (T cX (x)x{y)) ~ (x(x)) (x(y)) , (2.8) 

and reads 

T[G,D, X ] = S[0,x] + ^TrlnG- 1 + ^TrG^ 1 (G - G ) 






(2) 



(3) 



(4) 



Figure 2: NLO (2 loops) and N 2 LO (3 and 4 loops) contributions in the 2PI-l/iV 
expansion. The scalar propagator G is denoted with the full line and the auxiliary- 
field propagator D with the dashed line. 



+ l -TrlnD- 1 + ^Tr D l (D - D ) + F 2 [G,D}. (2.9) 

Since we take {<f> a ) = 0, there is no mixing between the and x propagators [THE]. 
The free inverse propagators read 



inverse propagators read 



3N 

G^ ab (x,y) = i[n x + m 2 + x(x)}5 ab 8 c (x-y), D \x,y) = —5 c (x-y). (2.10) 

The evolution equations for the propagators G and D and the one-point function \ 
are obtained by extremizing (|2.9|) and read 

_ p. + m , + m G-(I , B) = % J^ (X , j)0<l(z , y) + ^ _ B) , 

— D(s, y) = i J n(x, *)£>(*, y) + i<S c (a: - y), (2.11) 



and 



x(z) = gjy G «:(^a;). 



The self energies are defined by 

■5G„ t M5 

It is convenient to separate the local part of D J7j and write 



D(x,y) = -—- 
y ,yj 3N 



iS c (x-y) +D(x,y) 



such that D is determined from 
D(x,2/) = 



±Tl(x,y) + ^ fu(x,z)D(z,y). 



(2.12) 



(2.13) 



(2.14) 



(2.15) 
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(1) (2) (3) 

Figure 3: Self energy £ at NLO (1 loop) and N 2 LO (2 and 3 loops). 







(1) (2) (3) 

Figure 4: Auxiliary-field self energy II at NLO (1 loop) and N 2 LO (2 and 3 loops). 

We now continue with the 2PI-1/./V expansion, in which there is one diagram at 
NLO and two diagrams at N 2 LO, see fig. El For a detailed powercounting discussion 
we refer to ref . ^7] . It suffices here to say that a closed scalar propagator G ~ N and 
the auxiliary- field propagator D ~ 1/N. It is then easy to see that diagram (2) ~ 1 
and diagrams (3) and (4) ~ 1/N. The expressions are 

T^°[G,D] = - - I G 2 ab (x,y)D(x,y), (2.16) 

^ J xy 

[G,D] = -- / G ab {x, y)G bc (y, z)G cd {z, w)G da (w, x)D(x, z)D(y, w), 

" J XV zw 



^NNLO(3) 



^NNLO(4) [ 
2 



T" Ly ^>[G,D] 



Y2 



G ab (x, y)G bc (y, z)G ca (z, x) 



*G a , v tf, y')Gvc>{y', z')G c >a'(z', x')D(x, x')D(y, y')D(z, z'). 



The corresponding self energies are shown in figs. El and 0] 

We continue with the O(N) symmetric case, such that G a b(x,y) = G(x,y)5 a b and 
^ab{ x ,y) — ^(x, y)5ab. For notational simplicity we label the self energies according 
to the number of loops, e.g. S", where £ = 1,2,3. We stress that the 2PI-1/7V 
expansion does not coincide with loop expansion in the auxiliary-field formalism. For 
instance, there are two more four-loop diagrams, see fig. which only contribute 
at N 3 LO. The connection with the standard loop expansion is discussed further in 
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Figure 5: Loop expansion in the auxiliary- field formalism: additional four- loop dia- 
grams, of order 1/iV 2 (N 3 LO). 

appendix 1X1 

At NLO the self energies read 

E^(x,y) = -G(x,y)D(x,y), (2.17) 

n (1) (x,y) = -yG 2 (x,y), (2.18) 

and at N 2 LO 

rt 2) (x,y)= J G(x,w)G(w,z)G(z,y)D(x,z)D(w,y), (2.19) 

^ (x, y) = -N [ G(x, z)G(z, y)G(x', y')G(y', z')G{z' , x')D{x, x')D(y, y')D{z, z'), 



rt 2 \x,y) = j I G(x,w)G(w,y)G(x,z)G(z,y)D(z,w), ( 2 - 2 0) 



N 
I^\x,y) = / G(x,z)G(z,w)G(w,x)G(y,z')G(z',w')G(w',y)D(w,w')D(z,z'). 



3 Causal equations 

To bring the evolution equations in a form that can be solved numerically, the contour 
propagators and self energies are written in terms of statistical (F) and spectral (p) 
components [TUj . 

G{x, y) = F(x, y) - -sgn c (x, y)p(x, y), (3.1) 

and similar for D, E and II. Here sgn c (x, y) = Oe(x° — y°) — @c(y° — x°) is the sign- 
function along the contour, see appendix |Bj The explicitly causal equations then read 

na 

[D x + M 2 (x)] p{x, y) = - dz S p (x, z)p{z, y), (3.2) 

JyO 



r-X° py° 

[D x + M 2 (x)]F(x,y) = - dzY, p (x,z)F(z,y) + dzT, F (x,z)p{z,y), 

Jo Jo 



where 



M 2 (x) = m 2 + A — t^F(x, x), 
6N 



and 



3iV 

X 

3N 
X 



D P (x, y) = -Ii p (x, y)+ dz U p (x, z)D p (z, y), 

Jy° 

px° 

D F (x, y) = -U F (x, y)+ dz U p (x, z)D F (z, y) 



(3.3) 



(3.4) 



dz U F (x, z)D p (z, y). 



We use here the notation 



dz 



dz° / d 3 z. 



(3.5) 



Equations of motion derived from truncations of the 2PI effective action conserve the 
following energy functional (cf. eq. (J2.6JI ). 



E/N= J d 3 x- 



d x o d v o + d x i d v i + m 2 + 



±Z±l F {x,x) 



F(x,y) 



x=y 



+ - / d 3 x / dz° / d 3 z [T lp (x 1 z)F(z,x) —T lF (x, z)p(z,x)] . (3.6) 
The statistical and spectral self energies at NLO can be found in ref. [Ej and read 

E?>(*,y) = 



_A_ 

37V 
A 



F(x,y)D F {x,y) - -p(x,y)D p (x,y) 



3N 

4 1) (x,y) = -^ 



p(x, y)D F (x, y) + F(x, y)D p (x, y) 



F(x,y)F(x,y) - -p(x,y)p(x,y) 



U ( p 1 \x,y) = -NF{x,y)p{x,y). 



(3.7) 



We now come to the causal self energies at N 2 LO. These self energies have inter- 
nal vertices and therefore require further contour integrals. We start by discussing 



E( 2 )(x, y) in some detail. Since we have separated the local part of D, we first insert 
eq. (jTUll into eq. (12331) . This yields 

^ 2 \ X) y) = J D(x,z)G(x,w)G{y,z)D{y,w)G(z,w) 
= -g 2 G 3 (x,y) 



+ig< 



D(x, z)G(x, y)G 2 (y, z) + G 2 (x, z)G(x, y)D(y, z) 



where 



+g 2 / D(x,z)G(x,w)G(y,z)D(y,w)G(z,w), 



_\_ 

3iV' 



(3.8) 



(3.9) 



This is naturally organised according to the number of D propagators and we use the 
notation E^ 2 '™) for the contribution with n l)'s. The term without D propagators is 
the N 2 LO contribution to the setting-sun diagram and reads 



F&,y) = -g 2 



E 



F 2 (x,y)- -p 2 (x,y) 
3F 2 (x,y)- -p 2 (x,y) 



F(x,y), 
p(x,y). 



(3.10) 



The other terms quickly become rather lengthy, so we first discuss the general struc- 
ture. Every line in the self energy can be either a statistical (F) or a spectral (p) 
function. With L lines, this gives a maximum of 2 L possibilities. However, in order 
to find a nonzero result, every internal vertex needs at least one sgn- function coming 
from a p-type line ending on it (see appendix EI) • This implies that a diagram with 
V internal vertices should have at least V p-type lines. This reduces the maximal 
number of terms in the expressions for the causal diagrams to 



# distinct contributions < >^ ( , I < 2 J 



(3.11) 



where l p is the number of p-type lines. The actual number of nonzero contributions 
is in fact slightly less, since some contributions vanish after performing the contour 
integrals, due to the appearance of internal vertices without p-type lines (even though 
l p > V). Finally, the number of terms that have to be independently evaluated is 
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further reduced due to the fact that statistical (spectral) self energies are explicitly 
even (odd) under interchange of x and y. We also note that if the contour self energy 
is proportional to i (such as X^ 2,1 )), expressions with an odd number of p-type lines 
contribute to T, F and expressions with an even number of p-type lines to S p . If the 
contour self energy is proportional to 1 (such as E^ 2 ' 2 )), this is reversed. 

Using the contour integration rules summarized in appendix [0 we find explicitly 



-£ (2>1) ( 



x,y) 



+ dz 

Jo 
rv 
+ dz 

Jo 



2 F(.r :)p(x.z)D,,(,,. :)-[ DJx. :) ( F 2 {y, z) - -p 2 {y, z) 



2D F (x, z)F(y, z)p(y, z) + ( F 2 (x, z) - -p 2 (x, z)\ D p (y, z) 
D p (x, z)F(y, z)p(y, z) + F(x, z)p(x, z)D p (y, z) p(x, y), 



F(x,y) 
F(x,y) 

(3.12) 



and 



~ 2 ^\x,y) 

+ dz 

Jo 
rv 
+ dz 

Jo 



2F(x, z)p(x, z)D F (y, z) + D p (x, z) [ F 2 (y, z) - ^p 2 (y, z 



dz2 



p(x,y) 
p{x,y) 

D p (x,z)F(y,z)p(y,z) + F(x,z)p(x } z)D p (y } z) F(x,y), (3.13) 



2D F (x, z)F{y, z)p(y, z) + [ F 2 (x, z) - -p 2 {x, z) ) D p (y, z) 



and for the self energy with two D propagators 



£(2,2), 



■& F 



x,y) 



dz / dw p(x, z)D p (x,w)D F (y,z)F(y,w)F(z,w) 

ii Jo 

y ry 

dz / dw D F (x, z)F(x,w)p(y, z)D p (y,w)F(z,w) 
n Jo 

X PZ 

dz / dw p(x, z)D F {x 1 w)D F {y 1 z)F(y,w) 
n Jo L 

+D F (x, z)F{x, w)F(y, z)D p {y, w) p(z, w) 



y rz 
I dz dw 

<o Jo 



D F (x, z)F(x, w)p{y, z)D F (y, w) 
11 



+F(x, z)D p (x, w)D F (y, z)F(y, w) p(z, w) 



!l 



dz I dw 
o Jo 



D p (x, z)F(x, w)F(y, z)D p (y, w) 



+p(x, z)D F (x, w)D F (y, z)p(y, w 
y i 



F(z,w) 



dz I dwj D p (x,z)F(x,w)p(y,z)D p (y,w) 

p(x, z)D F (x, w)D p (y, z)p(y, w) p(z, w) 
D p (x, z)p(x, w)p(y, z)D F (y, w) 



y r-x | 

dz I dw - 
4L 



+p(x, z)D p (x, w)t)p(y, z)F{y, w) p(z, w) 



and 



+ dz dw-Dp(x,z)p(x,w)p(y,z)Dp(y,w)F(z,w), (3.14) 

Jy Jy 4 



T, p 2 ' 2) (x,y) = + dz dw D p (x,z)p(x,w)p(y,z)D F (y,w) 



+p(x, z)D p (x, w)Dp(y, z)F(y, w) 



F(z,w) 



+ I dz dw 

'y Jo 



D p (x, z)F(x, w)p(y, z)D p (y, w) 



+p(x, z)D F (x, w)Dp(y, z)p(y, w) 



F(z,w) 



+ dz dw 

ly Jo 



D p (x, z)F(x, w)p(y, z)D F (y, w) 



+p(x, z)D F (x, w)D p (y, z)F(y, w) p(z, w) 

X l>Z 1 

dz dw- D p (x,z)F(x,w)F(y,z)Dp(y,w) 
y Jy Zl 

+p(x, z)D F (x, w)D F (y, z)p(y, w) p(z, w) 

r x r z 1 r 
+ / dz dw- F(x,z)Dp(x,w)Dp(y,z)F(y,w) 

Jy J x *• '" 

+D F (x, z)p(x, w)p(y, z)D F (y, w) p(z, w) 
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x f y l r - 

dz dw- D p (x,z)p{x,w)p(y,z)D p (y,w) 
Jo 4 L 

+p(x,z)D p (x,w)D p (y,z)p(y,w) p(z,w). (3.15) 



In a few terms we used the symmetry of the integrand to make some minor simplifi- 
cations. 

For the auxiliary-field self energy we proceed in the same manner and find, at two 
loops, 



n (2) (*,</) 



N 



G(x, z)G(x, w)G(y, z)G(y, w)D(z, w) 



(3.16) 



' zw 

= l l [G 2 (x,z)G 2 (y,z) + ^ [ G(x,z)G(x,w)G(y,z)G(y,w)b(z,w). 
6 J z 6 J zw 

Denoting these diagrams again as n( 2,n ) , where n denotes the number of D propaga- 
tors, we find with n = 0, 

Ilf \x,y) = +± [ X dzF(x,z)p(x,z) (F 2 (y,z)-±p 2 (y,z) 



3 jo 
A rv 



1 



+ 3 / dz [ f2 ^ x,z " > ~~ ip 2 ( x > z n F (y> z )p(y> z )> 

2A f x 

Uf°\x, y) = -— J dz F{x, z)p(x, z)F(y, z)p(y, z), 



(3.17) 



and with one D propagator 



^F\x,v) 



"X rx 

dz I dw p(x, z)p(x,w)F(y, z)F(y,w)Dp(z,w) 
'o Jo 
ry rv 

— I dz dw F(x, z)F(x,w)p(y, z)p(y,w)D F (z,w) 
Jo Jo 

pz 

dz / dw2p(x, z)F(x,w)F(y, z)F(y,w)D p (z,w) 
'o Jo 

dz / dw2F(x, z)F(x,w)p(y, z)F(y,w)Dp(z,w) 
'o Jo 

— / dz j dw2p(x, z)F(x,w)F(y, z)p(y : w)Dp(z : w) 
Jo Jo 

rx rx -i 

+ I dz I dw -p(x, z)p(x,w)p(y, z)p(y : w)D F (z : w) 
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f y 1 
dz / dw -p(x, z)F(x, w)p(y, z)p(y, w)D p (z, w) 

[ x 1 
dz dw-p(x,z)p{x,w)p(y,z)F(y,w)D p (z,w), (3.18) 



and 



6 f x f x 

-Uf' 1] (x,y) = + dz dw2p(x,z)p(x,w)p{y,z)F(y,w)D F (z,w) 

r-x r-y 

+ l dz l dw2p(x, z)F(x,w)p(y, z)p(y,w)DF(z,w) 

Jy JO 



+ dz dw p(x, z)F(x,w)F(y, z)p(y,w)D p (z,w) 
■Jy -Jy 

PX PZ 

+ I dz dw F(x, z)p(x,w)p(y, z)F(y,w)D p (z,w) 



+ J dz dw2p(x, z)F(x,w)p(y, z)F(y,w)D p (z,w) 

rv l 
dz / dw-p(x,z)p(x,w)p(y,z)p(y,w)D p (z,w). (3.19) 
Jo * 

Since the II self energies have more internal symmetry than the £ self energies, the 
corresponding expressions are slightly shorter. We also emphasize that for nonequi- 
librium quantum fields statistical (F) and spectral (p) components are independent 
|1UJ : there are therefore no (obvious) cancellations between the various terms above. 

We now briefly discuss the requirements for a numerical solution. The evolution 
equations for the two-point functions (J3.2}) and (J3.4J) can be solved numerically by 
discretisation on a space-time lattice. For this it is necessary to perform the space- 
time integrals on the RHS of those equations. This has been discussed extensively in 
references cited in footnote El At NLO, the self energies (|3.7J) are simple products of 
F, p, Dp and D p . At N 2 LO however, the self energies themselves require space-time 
integrals to be performed: the two-loop self energy presented above contains up to 
two internal vertices, and the three-loop self energy will contain up to four internal 
vertices. These lead to nested loops over time, dramatically increasing the necessary 
CPU time. Because of the obvious numerical effort required, we have not attempted 
to solve the full N 2 LO approximation. Instead we concentrate in the next section on 
the three-loop approximation (two-loop self energy) in the + 1 dimensional case. 

The statistical and spectral self energies corresponding to the three-loop self en- 
ergies £( 3 ) and nJ 3 ) in eqs. (12.19)) and (J2.2U)) are very lengthy. For instance, for E^' 
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one has to consider 44 distinct nonzero diagrams which differ in the way the statistical 
and spectral functions appear on the 8 internal lines. Because we are not including 
these diagrams in the numerical analysis below, we refrain from giving the explicit 
expressions. Using the contour integrals listed in appendix [B] it is straightforward, 
albeit cumbersome, to work them out. 

4 Quantum dynamics 

Since a numerical solution of the full 3 + 1 dimensional field theory seems intractable, 
we now restrict our attention to quantum mechanics, which in this context is equiv- 
alent to 0+1 dimensional field theory. We consider the Hamiltonian 

ft = 2 PaPa + 2 m2qaqa + 4!iv ( <Ma ) 2 ' ( 41 ) 

and consider O(N) symmetric Gaussian initial conditions, parametrized by the nor- 
malized Gaussian density matrix [3H] . 



W\p(v,ti,v)\q)= (2<2)jV/2 e x P 



° + 1 (J2 , „2X , • V {,2 <l\ , V "I 



(la + C) + iirMa ~ %) + ^^T^a 



^£2 v±a iaJ 2£ 4£ 2 

(4-2) 

Recall that we consider (q a ) = (p a ) = 0. When a = 1, this reduces to the density 
matrix of a pure state, 

pfa,£,0) = |*o)<*o|, (4.3) 

with 

1 , .. V \ 2 



*o(g) = (g|*o) = (27re)jV /4 eX P 

In this case, the Schrodinger equation 



i2 +^ Qa 



4^ 2£ 



(4.4) 



i-y(q a ,t)=Hy(q a ,t), (4.5) 

can be solved numerically (see appendix |Hj), which allows for a comparison with the 
untruncated evolution |2T]. 5 

The expressions from the previous sections remain unchanged, provided that all 
reference to space indices and integrals are dropped, e.g., 

F ab ( Xj y) -> F ab (t,t') = 6 ab F(t,t') = -( qa (t)q b {t') + q b (t')q a (t)) , (4.6) 



3 In ref. |27| only density matrices with r/ = (and er = 1) were considered. 
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Pab(x,y) -» Pab{t,t') = S ab p(t,t') = i({q a (t),q b (t')}). 



(4.7) 



The density matrix ()4.2|) yields the following initial conditions for the statistical func- 
tion 



F(t,t') 
d t F(t,t') 

d t d t ,F(t,t') 



t=t'=0 



t=t'=o 






(U 



a 



t=t'=o 
As always, the spectral function satisfies 

p(t,t) = 0, d tP {t,t') 



"' + i?i- 



t=t> 



(4.9) 



The initial conditions for Dp and D p are determined by eq. ()3.4|) . The (conserved) 
energy takes the value 



(H)/N 



1 2 o^_ 
2 1 + %t? 



1 2 2 AiV + 2 4 

+ 2™* + 4!nv-^ 



(4.10) 



We solve the evolution equations (J3.2J) and (J3.4)) numerically for various A and AT 
using a simple leap-frog algorithm and a standard discretization of the time integrals. 
We set m = 1 throughout. For a step size of dt = 0.01 the entire memory kernel fits 
on a 1GB CPU, for runs until t = 20 (2000 timesteps). Using dt = 0.001 leads to 
indistinguishable results. A run which includes the first N 2 LO diagram takes about 2 
days on a single 3GHz machine, although presumably this can be improved somewhat. 
The CPU-time grows as the number of timesteps to the 4th power. A run at NLO is 
roughly 400 times faster. Including another nested integral (i.e. including the second 
N 2 LO diagram) is expected to give at least another factor of 400, making such an 
extension very challenging indeed. 

In order to study the effect of the different initial conditions, we show in fig. El 
the time evolution at NLO for N = 8 and A = 1, for various choices of initial states 
with the same energy. The early evolution, < t < 0.5, can be understood from the 
uncoupled case (A = 0), in which the dynamics is readily solved, and the two-point 
functions are 



F(t, t') = £ 2 cos(mt) cos(mt') 



1 



77T 



?/- 



a- 
4^ 



sin(mt) sin(mt') H sin [m(t + t )] , 



m 
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_ ^=0.4, r|=5.0, 0=1.0 

-- £=0.4, T)=4.5, 0=2.0 

£=0.2, ti=4.5, 0=1.0 

— £=0.4, ti=0.0, 0=4.1 
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Figure 6: Unequal time correlator F(t, 0) at NLO for various choices of initial states 
with the same energy (N = 8, A = 1). 



p(M') 



sin \m 



(t - t')\ 



m 



(4.11) 



The subsequent evolution is of course very different from the free case, in which 
damping is absent. 

We now continue with a relative large value of rj — 5 and a smaller £ = 0.4 
(corresponding to a squeezed initial state) since this yields an initially large amplitude 
and subsequent strong damping effects. For the sake of comparison with the numerical 
solution of the Schrodinger equation, we consider from now on pure states only, a = 1. 

In fig. [7|we show the evolution of the equal-time correlator F(t, t) for the different 
levels of truncation, at N = 8. In the Hartree approximation, all nonlocal terms on 
the RHS of the evolution equations (J3.2J) are dropped and only the time-dependent 
mass parameter is preserved. Eqs. (J3.4J1 are dropped altogether. At NLO the one- 
loop self energies (J3.7J) without internal vertices are kept. As explained at the end 
of the previous section, at N 2 LO we only keep the two-loop self energies, with up to 
two internal vertices, in addition to the one-loop self energies. We will refer to this 
truncation as the N 2 LO' approximation. Since it is derived from the two and three- 
loop diagrams in fig. |2l it is a 2PI self-consistent approximation. The result labelled 
with 'exact' corresponds to the numerical solution of the Schrodinger equation, which 
is detailed in appendix O Energy conservation at N 2 LO' is demonstrated in fig. |H1 
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Figure 7: Time evolution of the equal-time correlator F(t, t) for various approxi- 
mations and the exact result from a numerical solution of the Schrodinger equation 



15r 



— total 

— kinetic 
potential 

— memory term 




Figure 8: Time evolution of different energy components at N 2 LO' (same parameters 
as in fig. Ej) . 
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where the three components in eq. (|3.6|) and their sum is shown. 

A close look at fig. UJ shows that for early times both the NLO and the N 2 LO' 
approximation are in quantitative agreement with the exact result. The amplitude in 
the Hartree approximation shows no sign of decreasing, due to a complete absence of 
dephasing in quantum mechanics. 6 Around t = 2.7, the NLO approximation starts to 
differ from the exact evolution, whereas the evolution at N 2 LO' is capable to follow 
the exact evolution a bit longer. Around t = 4 we find that both truncations fail 
to track the exact evolution and continue to evolve in an irregular fashion. We have 
verified that this behaviour is not due to the time discretisation. We also note that 
energy remains conserved. The irregular behaviour at later times seems to be peculiar 
to quantum mechanics and has been observed before in dynamics from truncated 
effective actions [27] . As far as we know it has not been observed in 2PI dynamics in 
the field theory case, where already the NLO approximation results in equilibration 
and thermalization. In this sense quantum mechanics, with only a finite number of 
degrees of freedom, is very different. 

If we continue with a comparison at early times, we note that for N = 8 the 
N 2 LO' approximation works slightly better than the NLO approximation. This can 
be investigated further by looking at different values for N. Since the N 2 LO' term 
is suppressed by 1/N, it is expected that the difference between the N 2 LO' and the 
NLO evolution will be largest for smaller N, while for larger N the 1/N expansion 
itself is better behaved and the difference between N 2 LO' and NLO is reduced. 

This qualitative picture is confirmed by first going to smaller N. In fig. E3 we 
show again the equal-time correlation function, but now for N = 2. As expected, the 
evolution ceases to follow the exact one earlier, but a close look at the first maximum 
indicates that it is first the Hartree approximation that breaks down, subsequently the 
NLO approximation and finally the N 2 LO' approximation. At later times irregular 
behaviour is again observed. We find therefore that at early times an increase in the 
order of the truncation has a quantitatively correct effect. In fig. [TOlthe equal-time 
correlation function is shown again, but now for larger N = 20. In this case the effect 
of adding the N 2 LO' contribution is much less important, as expected, and we find 
that the N 2 LO' evolution follows NLO rather than the exact curve, consistent with 
an effective convergence of the expansion for large N. Both curves follow the exact 
evolution for longer than in the cases shown above. 

A comparison between the NLO, the N 2 LO' and the exact evolution is shown in 



6 This is different in field theory, where the amplitudes of equal-time correlation functions diminish 
in time due to dephasing. However, in the Hartree approximation the system still experiences 
recurrence. This shows that damping in the Gaussian approximation is not a result of irreversible 
loss of memory. 
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— Hartree 
-- NLO 

— NNLO' 
exact 




Figure 9: Same as Figured for N = 2. 



— Hartree 
-- NLO 

— NNLO' 
exact 




8 10 

Figure 10: Same as Figured for N = 20. 
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Figure 11: Time evolution of F(t,t) for NLO (upper), N 2 LO' (middle) and exact 
(lower) dynamics for various 2 < iV < 20 (other parameters as in fig. d)- 
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Figure 12: Unequal-time correlation function F(t, 0) for NLO (upper) and N 2 LO' 
(lower) dynamics for various 2 < N < 20 (other parameters as in fig. Ej). 
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figure ITT1 for a wide range of N, from A" = 2 to iV = 20. It is observed that the exact 
solution appears to be intermediate between the NLO and N 2 LO' evolution, with the 
evolution at N 2 LO' performing slightly better, in particular in terms of the oscillation 
frequency at smaller N. 

As mentioned above, thermalization and equilibration cannot be investigated in 
quantum mechanics. However, in analogy with field theory effective loss of memory 
can be studied. In fig. ED we show the unequal-time two-point function F(t, 0) for 
various values of N. Because it is nonlocal in time, this correlator is not immediately 
accessible from the solution of the Schrodinger equation. Smaller values of A" corre- 
spond to a more rapid decrease of the amplitude, as expected. Maybe surprisingly, 
the memory appears to be washed out faster at NLO than at N 2 LO' . However, this 
conclusion should be treated with care, since the approximations fail to track the ex- 
act evolution after some (N dependent) time, as was shown above for the equal-time 
correlation functions. 

5 Outlook 

We considered nonequilibrium dynamics in the O(N) model, employing the 2PI-1/AT 
expansion to N 2 LO. We presented the explicit expressions for the three- loop approx- 
imation in the auxiliary-field formalism, to which we refer as N 2 LO', and indicated 
how to obtain to full N 2 LO contribution. The resulting evolution equations were 
solved numerically for quantum mechanics. 

While the qualitative change in the nonequilibrium evolution when going from 
Gaussian (LO) approximations to NLO is enormous, the impact when changing from 
NLO to N 2 LO is reassuringly small. This indicates that the 2PI-1/A" expansion is ef- 
fectively rapidly converging and that higher order effects give quantitative corrections 
only. We found that at early times the evolution at N 2 LO' performs slightly better 
than at NLO. This is especially visible at small N, where the higher-order contribu- 
tion is not suppressed. At large times we found that all truncations break down, but 
we believe that this is special for quantum mechanics since it has not been observed 
in field theory, where already the NLO approximation has been seen to perform very 
well. It would therefore be very interesting to implement the N 2 LO truncation, or 
at least the three-loop diagram as we considered here, in a 1 + 1 dimensional scalar 
O(N) model and test the apparent convergence quantitatively in field theory, also for 
late times. Alternatively, with some effort one may also be able to do the full N 2 LO 
approximation in the case of quantum mechanics. 

Finally, in this paper we only considered the symmetric phase. However, an 
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extension to the broken phase at N 2 LO is straightforward, since it involves only one 
additional two-loop diagram in the effective action, yielding new self energy type 
contributions without internal vertices ^7]. In fact, this diagram has been included 
already in the so-called 'bare vertex approximation' (BVA) [Tlj . 

Acknowledgments. G.A. is supported by a PPARC Advanced Fellowship. A.T. is 
supported by the PPARC SPG "Classical Lattice Field Theory" . 



A Loop expansion 

In order to make a connection between the 1/iV expansion and the ordinary loop 
expansion (see fig. IT5)I . we give here the expression up to five loops in the loop 
expansion in the 0(A)-symmetric case. We write the 2PI contribution as T2 = : 

e: 



^ 2 T 2 and find to fourth order 



,(2) 



,(3) 



,("1) 



A (A + 2) 

8 3 
i\ 2 (N + 2) 
48" 3A 



G (x, x) 



G 4 (x,y), 



XII 



A 3 (A + 2)(A + 8) 



48 



27A 2 



G 2 (x,y)G 2 (x,z)G 2 (z,y). 



xyz 



At fifth order two diagrams contribute, T 2 



(5) 



,(5a) 



+ r 



(56) 
2 5 



which read 



r 



x 



(5a) ^ (A + 2)(A 2 + 6A + 20) 

128 81A 3 

G 2 {x, y)G 2 {y, z)G 2 {z, w)G 2 {w, x) 

xyzw 

, (56) a 4 (A + 2)(5A + 22) 
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S1A 3 



(A.l) 
(A.2) 
(A.3) 



(A.4) 




Figure 13: Loop expansion: contribution to T 2 up to fifth loop order. 
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Figure 14: Keldysh contour in the complex-time plane. 



G 2 {x, y)G(x, z)G{x, w)G 2 (z, w)G(y, z)G{y, w) 



(A.5) 



xyzw 



<2) 



In the 2PI-1/N expansion to N 2 LO, T 2 and T 2 are completely included, while from 



(3) 



<*) 



i(5o) 



(56) 



ri, and Tg the NLO and N LO parts are taken into account. The eye diagram T 2 ' 
starts at N 2 LO only. Part of the leading N 2 LO contribution is taken into account 
via the three-loop diagram in the auxiliary-field formalism and part via the four-loop 
diagram. The eye diagram is special since it is the first diagram that contributes to 
the bulk viscosity in the weak coupling limit 



B Contour integrals 

The evolution equations are formulated along the Keldysh contour in the complex- 
time plane, see fig. El Splitting the propagators and self energies in statistical and 
spectral components yields the real and causal equations discussed in sectional While 
for the NLO approximation this procedure is straightforward, it becomes cumbersome 
when multi-loop contour integrals are encountered and integrals over products of sgn- 
functions have to be evaluated. Here we give some general expressions we find useful. 
First consider integrals over products of B-functions. Since the contribution from 
the upper and lower part of the contour differ only in sign, one finds e.g. that 



[dz° = 0, I dz°e c (x,z)= [ dz°, 

Jc Jc Jo 



(b.i; 



or in general 



„ M M 

dz°l[e c (xlz°) = J2 






dz°. 



(B.2) 



Here we have taken the initial time at z° = 0. These results can be employed for 
integrals over products of sgn- functions, using that sgn c (x, y) = 2Q C (x, y) — 1, yielding 
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e.g. 



[ dz°sgn c {x,z) = 2 [ dz°, (B.3) 

Jc ' Jo 

or in general 



M M 



[dz°f[sgn c (xlz )=2<r Hsgndx^xl) H dz°. (B.4) 

In a theory with quartic interactions, ()B.4|) is needed with M < 4, for which we have 
verified this identity explicitly. 

C Numerical solution of the Schrodinger equation 

In the case of quantum mechanics, we have the option to solve the Schrodinger equa- 
tion numerically for the full wavefunction. This only applies to pure states, a — 1. 
The system corresponds to a spherically symmetric anharmonic oscillator in N di- 
mensions, and by imposing O(N) symmetry and make suitable redefinitions we can 
reduce the problem to that of a particle on a halfline in one dimension 27J. 

If the wave function is written as the product of a radial function, depending on the 
(rescaled) radial coordinate < r = \/q a q a /N < oo, and hyperspherical harmonics, 
depending on the N — 1 angles, 

*(g a ,t) = (Nr 2 ) {N - 1)/4 ^(r,t)Y(Q), (C.l) 

the O(N) symmetric problem is reduced to the effective one- dimensional Schrodinger 
equation [27] 

^*{r,t)=n*$(r,t), (c.2) 

with the Hamiltonian 



^ = -mw +uir) > (c - 3) 



and the effective potential 



„ M (l-l/iV)(l-3/iV) 1 22 A 4 ,_,, 

U i r ) = ~ —& — + 2 mr + 4! r ■ (C ' 4) 
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For a pure state the initial radial wave function, corresponding to the density matrix 
discussed in section HJ is 



i / N\ N/4 

<I>(,-.0): (^L) r^)/ 2 exp 



V \ _2 



-JV — + z — r 
A? 2^ 



(C.5) 



This initial wave function is normalized with respect to the innerproduct 

POO 

($|$) = / dr$\r,t)$(r,t). (C.6) 

Jo 



We solve eq. (JC.2J) numerically using the second-order Crank-Nicholson differencing 
scheme J31]. For the time intervals shown, both unitarity and energy are preserved 
better than 1 in 10 12 . 



References 

[1] For a recent review, see J. Berges, AIP Conf. Proc. 739 (2005) 3 
|hep-ph/0409233 . 



[2] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena, Oxford Univer- 
sity Press (1989); H. J. de Vega, Commun. Math. Phys. 70 (1979) 29; J. Avan 
and H. J. de Vega, Phys. Rev. D 29 (1984) 2904. 

[3] For examples of nonequilibrium quantum field dynamics in the mean-field ap- 
proximation, see e.g., F. Cooper, S. Habib, Y. Kluger, E. Mottola, J. P. Paz 
and P. R. Anderson, Phys. Rev. D 50 (1994) 2848 [hep-ph/940 5352|; D. Boy- 
anovsky, H. J. de Vega, R. Holman, D. S. Lee and A. Singh, Phys. Rev. D 51 
(1995) 4419 |hep-ph /9408214|; D. Boyanovsk y, H. J. de Vega, R. Holman and 
J. F. Salgado, Phys. Rev. D 54 (1996) 7570 |he p-ph/9608205] ; D. Boyanovsky, 
C. D estri, H. J. de Ve ga, R. Holman and J. Salgado, Phys. Rev. D 57 (1998) 
7388 |hep-ph/9711384|; J. Baa cke, K. Heitmann and C. Patzold, Phys. Rev. D 
58 (1998) 125013 |hep-ph/ 9806205|, and references therein. 

[4] L. M. A. Bette ncourt and C. Wetterich, Phys. Lett. B 430 (1998) 140 
|hep-ph/9712429| . 

[5] G. Aarts, G. F . Bonini and C. Wetterich, Phys. Rev. D 63 (2001) 025012 
|hep-ph/0007357 . 



27 



[6] G. Aarts and J. Smit, Phys. Rev. D 61 (2000) 025002 |hep^ h/9906538] . 

[7] M. Salle, J. Sm it and J. C. Vink, Phys. Rev. D 64 (2001) 025016 
|he p-ph/0012346| . 

L. M. A. Bettencou rt, K. Pao and J. G. Sanderson, Phys. Rev. D 65 (2002) 
025015 hep-ph/0104210|. 



[9] J. Berges and J. Cox, Phys. Lett. B 517 (2001) 369 |hep-ph/0006160|. 

[10] G. Aarts and J. Berges, Phys. Rev. D 64 (2001) 105010 lhep-ph/0103049"l . 

[11] S. Juchem, W . Cassing and C. Greiner, Ph ys. Rev. D 69 (2004) 025006 
|he p-ph/03073 53|; Nucl. Phys. A 743 (2004) 92 |nucl-th/0401046| . 

[12] A. Arrizaba laga, J. Smit and A. Tranberg, Phys. Rev. D 72 (2005) 025014 
|hep-ph/0503287|. 



[13] J. Berges, Nucl. Phys. A 699 (2002) 847 hep-ph/0105311 



[14] B. Mihaila, F. Cooper and J. F. Dawson, Phys. Rev. D 63 (2001) 096003 
|hep-ph/0006254| . 

[15] G. Aarts and J. Berges, Phys. Rev. Lett. 88 (2002) 041603 |hep-ph/0107l"29| . 

[16] K. Blag oev, F. Co oper, J. Dawson and B. Mihaila, Phys. Rev. D 64 (2001) 
125003 |hep-ph/0106l95l; F. Cooper, J. F. Dawson and B. Mihaila, Phys. Rev. 
D 67 (2003) 051901 |hep-ph/020"7346"l . 

[17] G. Aarts, D. Ahrensmeier, R. Baier, J. Berges and J. Serreau, Phys. Rev. D 66 
(2002) 045008 |hep-ph/0201308 . 



[18] J. Berges and J. Serreau, Phys. Rev. Lett. 91 (2003) 111601 hep-ph/0208070 



[19] F. Cooper, J. F . Dawson and B. Mihaila, Phys. Rev. D 67 (2003) 056003 
|hep-ph/0209051| . 

[20] A. Arrizab alaga, J. Smit and A. Tranberg, JHEP 0410 (2004) 017 
|hep-ph/0409177| . 

[21] J. Berges, S. B orsanyi and J. Serreau, Nucl. Phys. B 660 (2003) 51 
hep-ph/0212404|; J. Berges, S. Borsanyi and C. Wetterich, Phys. Rev. Lett. 
93 (2004) 142002 |hep-ph/0403234| . 

28 



[22] E. Calzetta and B. L. Hu, Phys. Rev. D 37 (1988) 2878. 
[23] M. Lindner and M. M. Muller, |hep-ph/0512147| 



[24] J. Berges and S. Borsanyi, hep-ph/0512155 



[25] E. A. Calzetta, B . L. Hu and S. A. Ramsey, Phys. Rev. D 61 (2000) 125013 
|hep-ph/9910334| . 

[26] G. Aarts and J. M. Martinez Resco, Phys . Rev. D 68 (2003) 085009 
|he p-ph/03032 16|; JHEP 0402 (2004) 061 |h ep-ph/0402192| ; JHEP 0503 (2005) 
074 |hep-ph/050316J] . 

[27] B. Mihaila, T. Athan, F. Cooper, J. Dawson and S. Habib, Phys. Rev. D 62 
(2000) 125015 jhep-ph/0003T05|; B Mihaila, J. F. Dawson, F. Cooper, M. Brew- 
ster and S. Habib, |hep-ph/9808234. 



[28] J. Berges and I. O. Stamatescu, Phys. Rev. Lett. 95 (2005) 202003 
|hep-lat705 08030|. 

[29] J. M. Cornwall, R. Jackiw and E. Tomboulis, Phys. Rev. D 10 (1974) 2428. 

[30] S. R. Coleman, R. Jackiw and H. D. Politzer, Phys. Rev. D 10 (1974) 2491. 

[31] L. F. Abbott, J. S. Kang and H. J. Schnitzer, Phys. Rev. D 13 (1976) 2212. 



[32] S. Jeon, Phys. Rev. D 52 (199 5) 3591 [hep -ph/9409250|; S. Jeon and L. G. Yaffe, 
Phys. Rev. D 53 (1996) 5799 |hep-ph/9512263| . 

[33] F. Cooper, S. Ha bib, Y. Kluger and E. Mottola, Phys. Rev. D 55 (1997) 6471 
|he p^ph/9610345| . 

[34] W.H. Press, S.A. Teukolsky, W.T. Vetterling and B.P. Flannery, Numerical 
Recipes in C, Second Edition (Cambridge University Press, United Kingdom). 



29 



